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Abstract. The polarization curve of GRB 020813 is discussed and compared to different models for the structure, evolution 
and magnetisation properties of the jet and the interstellar medium onto which the fireball impacts. GRB 020813 is best suited 
for this kind of analysis for the smoothness of its afterglow light curve, ensuring the applicability of current models. The 
polarization dataset allows us to rule out the standard GRB jet, in which the energy and Lorentz factor have a well defined 
value inside the jet opening angle and the magnetic field is generated at the shock front. We explore alternative models finding 
that a structured jet or a jet with a toroidal component of the magnetic field can fit equally well the polarization curve. Stronger 
conclusions cannot be drawn due to the incomplete sampling of the polarization curve. A more dense sampling, especially at 
early times, is required to pin down the structure of the jet and the geometry of its magnetic field. 



Key words, gamma rays: bursts - polarization - radiation mechanisms: non-thermal 



1. Introduction 



Send offprint requests to: D. Lazzati, lazzatiSast. cam . ac . uk 
* Based on observations collected at the European Southern 
Observatory, Cerro Paranal (Chile), ESO programmes 69.D-0461(A) 
and69.D-0701(A). 



The discovery that Gamma-Ray Burst (hereafter GRB) after- 
glows are linearly polarized is one of the strongest proofs that 
the photons we observe are synchrotron radiation (Covino et 
al. 1999; Wijers et al. 1999). Despite that, polarimetric after- 
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glow observations have so far been sparse and discontinuous, 
with only few polarimetric measurements performed on differ- 
ent optical transients (see Covino et al. 2003b for a review). 
Recently, however, several better sampled linear polarization 
measurements allowed for the first studies of the evolution 
of polarization, especially for GRB 021004 (Rol et al. 2003; 
Lazzati et al. 2003) and GRB 030329 (Greiner et al. 2003). 

These studies (Lazzati et al. 2003; Nakar & Oren 2003), 
coupled with theoretical works (Ghisellini & Lazzati 1999; 
Sari 1999; Rossi et al. 2004; Granot & Konigl 2003) reached 
two important conclusions. First, polarimetric studies provide 
unique information on the structure and dynamics of GRB out- 
flows: the polarization, e.g., from a homogeneous jet has a 
markedly different evolution from that of a structured one, even 
though their light-curves are barely distinguishable (Rossi et al. 
2004). Second, polarization studies are complex and subject to 
systematic errors. This is due to the combination of the intrin- 
sic weakness of GRB polarization ( <, 3%, at least at the times 
at which measurements have so far been possible; Covino et al. 
2003b, but see also Bersier et al. 2003 who found a possible 
polarization flickering at the 10% level) and to the sensitivity 
of the polarization signal to bright spots on the fireball surface 
and/or inhomogeneities in the ambient medium (Lazzati et al. 
2003; Granot & Konigl 2003; Nakar & Oren 2003). In addition, 
the measured polarized signal is comparable to the polarization 
induced by the propagation of light through a moderately ab- 
sorbing interstellar material, making the direct comparison of 
models with data more difficult (Lazzati et al. 2003). 

All these considerations make the dataset obtained for 
GRB 020813 unique and extremely interesting. GRB 020813 
(see the accompanying paper, Gorosabel et al. 2003, for more 
details) had an extremely smooth light curve (Laursen & 
Stanek 2003; Gorosabel et al. 2003), successfully fitted by a 
smoothly broken power-law with an r.m.s. scatter of the residu- 
als of < 0.01 magnitudes in the optical filters. This ensures that 
inhomogeneities in the fireball structure and/or in the surround- 
ing interstellar medium (ISM) are not significant and therefore 
can not affect the polarization measurement. In addition, the 
spectropolarimetric measurement of Barth et al. (2003) does 
not show evidence of a strong colour dependence of the polar- 
ization, a signature of the polarization induced by the interstel- 
lar medium (Serkowski et al. 1975). 

In this paper we compare the polarization evolution of 
GRB 020813 with existing models from the production of po- 
larized light in spherical and beamed fireballs. Some of these 
models, which were originally computed in a uniform environ- 
ment, are extended to the wind case (Sect. 2). We also consider 
the possible presence of an ordered component of the magnetic 
field advected from the central source (Sect. 2). The compari- 
son of the models with the data is described in Sect. 3 and we 
finally discuss our findings in Sect. 4. 

2. The models 

In this section we describe the models we will compare to the 
polarimetric data of GRB 020813. Although some of the theo- 
retical polarization curves are taken from the literature those as- 
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Fig. 1. Polarization curves from a homogeneous non sideways 
expanding jet with shock generated magnetic field in a stan- 
dard ISM and wind environment. Both curves are computed 
for observers located at 6> Q b s = 0.676^. The break time tj e t 
has been measured independently for the two light curves. The 
polarization position angle rotates by 90° in the moment of null 
polarization. 



sociated with magnetised jets are originally calculated for this 
work. 

First, let us consider the magnetic domain model, the first 
ever considered model for the observation of linear polariza- 
tion in GRB afterglows (Gruzinov & Waxman 1999). If the 
shock generated field is able to rearrange rapidly in ordered 
domains, Gruzinov & Waxman (1999) calculated that an aver- 
age observer should see about ~ 50 magnetic domains, and 
therefore if each domain produces a polarization po a net po- 
larization n = po/y/N ~ 0.1(p /0.7)(50/A^) 1 / 2 should be 
observed, where po ~ 70% for synchrotron. Deriving a po- 
larization curve for this model is not an easy task, given its 
intrinsic random character. Nevertheless, a general conclusion 
can be drawn. Polarization should be variable and variability in 
the degree of polarization should be associated to variability in 
the position angle (Gruzinov & Waxman 1999). In particular, 
a variation of linear polarization by a factor of 2 should be as- 
sociated to a random re-shuffling of the position angle. This is 
not what we observe in the data, where only a moderate rota- 
tion of a few degrees is associated to the polarization evolution 
(Gorosabel et al. 2003). We conclude therefore, analogously to 
what derived by Barth et al. (2003) (see also Greiner at al. 2003 
for the case of GRB 030329), that random patches of ordered 
magnetic field are not producing the observed polarization. 
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Fig. 2. Polarization curves for an MHJ jet with toroidal 
magnetic field. From left to right different curves refer to 
(Us/fljct = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9. The po- 
sition angle of the polarization is constant throughout the entire 
evolution. 



2. 1. Hydrodynamic Homogeneous Jet 

We define a "Hydrodynamic Homogeneous Jet" (HHJ) as a 
standard top-hat jet in which the energy, which is uniformly 
distributed within the jet opening angle, is carried by baryons 
and the magnetic field responsible for the afterglow syn- 
chrotron emission is tangled on small scales but overall domi- 
nated by either a component orthogonal or parallel to the shock 
front. Such a field configuration is generated either by compres- 
sion of a fully tangled magnetic field of by two stream plasma 
instabilities at the shock front, for example the Weibel instabil- 
ity (Silva et al. 2003). Polarization curves from this class of jets 
have been computed by various authors (Ghisellini & Lazzati 
1999; Sari 1999; Granot & Konigl 2003; Salmonson 2003; 
Rossi et al. 2004). Most of these papers, and in particular the 
more recent ones based on numerical computations rather than 
on analytical approximations, agree qualitatively on the result- 
ing polarization curve: it has two peaks, separated by a moment 
of null polarization roughly coincident with the break time of 
the total light curve. In this moment the position angle of the 
polarization rotates by 90°. The second peak is always stronger 
than the first, their ratio depending on the dynamics of the jet 
sideways expansion (SE): the faster the expansion, the smaller 
the second peak (the first peak is obviously only marginally 
affected by sideways expansion; see Rossi et al. 2004). In this 
paper we adopt the polarization curves computed by Rossi et al. 
(2004). After the jet break time we consider either no lateral ex- 
pansion or a jet expanding at the speed of sound in the shocked 
fluid comoving frame (Eq. 8 of Rossi et al. 2004). These two 
assumptions are chosen in order to encompass the numerical 
results of Kumar & Granot (2003), who find a sub-sonic ex- 



pansion until the expansion velocity becomes trans-relativistic. 
Extremely high sideways expanding efficiency were instead as- 
sumed by Sari (1999). In this case a different behaviour of the 
post-break polarization is obtained, where the polarization can 
have either no change in position angle or a double change (see 
also Barth et al. 2003). The energy density within the jet open- 
ing angle is assumed to be uniform before and after the break 
time. 

In addition to the calculations presented in the above men- 
tioned papers, we show here the effect of a wind environment 
on the polarization curve. In Fig. [2 a polarization curve for a 
HHJ expanding in an ISM environment (solid line) is compared 
to that of the very same jet that propagates in a wind (where 
n(r) oc r -2 ; dashed line). The qualitative result is identical, 
with the polarization curve characterised by two peaks with 
position angle shifted by 90°. The only difference is that in the 
wind case the evolution is slower, as already noted by Kumar & 
Panaitescu (2000) for the break of the total light curve. Before 
the non-relativistic transition (t <C ^nr) it is possible to obtain 
the wind polarization curve with good approximation from the 

3/2 

ISM one rescaling the times according to < w ; n d = ^ism- This, 
for a non sideways expanding jet, comes from the fact that the 
polarization behaviour mainly depends on the Lorentz factor T 
which scales as f ~ 3 / 8 for the ISM case and as < -1 / 4 in the wind 
environment. 

2.2. Hydrodynamic Structured Jet 

A "Hydrodynamic Structured Jet" (HSJ) is similar to a HHJ 
from the micro-physical point of view, but has a distribution 
of energy per unit solid angle which is larger in the centre and 
smaller in the wings (Rossi et al. 2002). In order to preserve 
the standard jet energy (Frail et al. 2001), this distribution must 
have the form of a power-law with the energy scaling as 8~ 2 . 
Polarization curves for this jet have been computed by Rossi 
et al. (2002; 2004). They are single peaked, with the time of 
the maximum coincident with the break time of the total light 
curve. The position angle does not vary throughout the evolu- 
tion, and the maximum observed polarization grows with the 
angle that the line of sight makes with the jet axis. For this 
case we did not consider sideways expansion since it affects 
only marginally the behaviour of the polarization curve, to a 
level much smaller than the accuracy of the data we are dealing 
with. A discussion of theoretical light curves can be found in 
Rossi et al. (2004). 

2.3. Magnetised Homogeneous Jet 

We call here "Magnetised Homogeneous Jet" (MHJ) a jet in 
which the magnetic field has a toroidal structure, with the po- 
lar axis coincident with the jet axis. The energy is uniform 
within the jet opening angle, analogously to a HHJ. Such a 
magnetic configuration can be realized, e.g., if the energy is 
carried in rough equipartition by the baryons and an electro- 
magnetic component (e.g. Proga et al. 2003). Such a scenario 
is naturally envisaged in the original GRB jet, where a siz- 
able magnetic field component can be advected from the cen- 
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Fig. 3. Same as Fig. |2]but with the time shown in units of 
tobs (defined through r(< Q b s ) = 9^ s ) rat her than t ]e t (defined 
throughr(t jet )=0r t 1 ). 

tral engine. Whether this field can be propagated to the external 
shock is less certain, but such a possibility should be considered 
(Lyutikov, Pariev & Blandford 2003; Lyutikov & Blandford 
2004), given the fact that contact discontinuities are known to 
be maximally unstable (Landau & Lifshitz 1989) due to the 
lack of any restoring force. Here we assume that all the syn- 
chrotron emission is produced by the toroidal magnetic field, 
ignoring a possible (even likely, given the low measured levels 
of polarization) random component of the field (see below). We 
also consider uniform jets, i.e. jets with a well-defined cone an- 
gle within which the jet is uniform. Structured magnetised jets 
will be discussed in the next section. Polarization curves are 
shown in Figs.[2]and[3] 

There are clearly major differences with respect to the hy- 
drodynamic jets (either homogeneous or not). The first is that 
the polarization does not disappear for very small observer 
times. This is due to the fact that in the hydrodynamic jets the 
polarization is cancelled by symmetry unless the edge of the 
jet (for the HHJ) or the brighter core (HSJ) are visible. In the 
MHJ case, the observed afterglow polarization is not due to a 
lack of cancellation, but rather to the presence of a genuinely 
ordered magnetic field. The higher the Lorentz factor of the 
fireball (and therefore the smaller the observer time) the less 
curved is the observed field on the plane of the sky and there- 
fore the higher the polarization (see e.g. Lyutikov et al. 2003). 
For ^obs/ #jet < 0.6 the polarization curve is monotonically de- 
creasing, while for 9 b s /9 iC t > 0.6 there is a maximum in the 
polarization curve at t ~ £ b s , where i b s is defined through 
r(iobs) = l/#obs- The position angle is constant throughout 
the whole evolution. 

Another important difference between hydrodynamic po- 
larization curves and MHJ ones is that in the case of hydrody- 
namic light and polarization-curves there is only one relevant 



time-scale, while MHJ curves have two time-scales. As can be 
seen by comparing Fig. [2] with Fig.|5J the break in the polar- 
ization curve takes place approximately at t = i b s , while the 
break in the total light curve happens at t = tj Ct . In the case 
of hydrodynamic jets, the relevant change of behaviour in the 
polarization curve (the change of angle in HHJs or the peak in 
HS Js) takes place in coincidence with the change of slope in the 
total light curve, i.e. the break time tj Ct . As a consequence, in an 
MHJ jet the simultaneous observation of light and polarization- 
curves can completely solve the jet and the observed geometry, 
something that is impossible in a hydrodynamic jet. 

Before discussing the last configuration, that of a force-free 
magnetic bubble, it is worth stressing that the curves shown in 
Figs.|2]and[3]heavily rely on the assumption that the magnetic 
field contained in the ejecta is toroidal. Such a field should be 
transported out to the shocked ISM producing the afterglow 
radiation. In our computations it is assumed that all the mag- 
netic field present in the shocked ISM comes from the GRB 
ejecta. It is more likely that a turbulent magnetic field is gen- 
erated at the shock and mixed with the ordered magnetic field 
of the ejecta. This is also required by the fact that the typi- 
cal afterglow polarization is at the level of few per cent, much 
smaller than the values predicted in Figs. [2] and [3] The sim- 
plest assumption that can be made is that the ratio between the 
turbulent and ordered components of the field stays constant 
throughout the evolution. In this case the resulting polariza- 
tion curve may be obtained by rescaling those in the figures. 
However, the ratio between the two fields may change with 
time and in that case, even if it is likely that the general mono- 
tonic behaviour would be maintained, the actual shape of the 
curves may change. Given the quality of the dataset and the 
theoretical difficulties in the propagation of the field, we adopt 
a completely ordered field or a mix with constant ratio. 

2.4. Magnetised Structured Jet 

It has been proposed recently that many astrophysical jets may 
be dominated by electromagnetic forms of energy rather than 
by baryonic matter (Lyutikov & Blandford 2004; Lyutikov et 
al. 2003). If the plasma is sufficiently tenuous, the flow can be 
followed under the Force-Free approximation, a set of equa- 
tions in which the inertia of the matter is neglected. We here 
consider the late phase of the evolution, when the external 
shock has developed and its dynamics does not differ any more 
from that of an hydrodynamic jet. Analogously to the MHJ jet, 
we assume that the magnetic field in the shocked ISM is trans- 
ported from the magnetic bubble. We call this configuration a 
"Magnetised Structured Jet" (MSJ), i.e. a jet in which the en- 
ergy distribution is inhomogeneous as described for the HSJ, 
and the magnetic field is toroidal as in the case of a MHJ. The 
difference with respect to the MHJ jet is therefore that the jet 
is not uniform within its opening angle but structured, with an 
energy per unit solid angle distribution Eq oc 1/ sin 2 9 ~ 9~ 2 , 
analogously to the HSJ. The polarization curve has been com- 
puted From Eq. 2, 6 and 7 of Lyutikov et al. (2003) and is 
shown in Fig.|4] The polarization curve for this class of jets is 
unique, independent of the observer line of sight 6> bs, if plot- 
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Fig. 4. Polarization curve for a MS J with purely toroidal mag- 
netic field. The polarization curve is unique and independent of 
the observer line of sight as long as the time is plotted in unit 
of the light curve break time. The position angle of the polar- 
ization is constant throughout the entire evolution. 

ted against the observer time in units of the jet break time tj Ct . 
The position angle is independent of time, since reflects the 
orientation of the magnetic field. As for the MHJ jet, a random 
magnetic field component is likely to be mixed in the toroidal 
field altering the detailed shape of the light curve, but not its 
general behaviour. 

3. Comparison of the data with the models 

The polarization dataset for GRB 020813 has been presented 
elsewhere (Gorosabel et al. 2003) and we here simply note two 
differences between their data and the data we adopt. First, we 
adopt the dataset uncorrected for polarization induced by the 
propagation into the Galactic ISM, as derived from the polar- 
ization properties of field stars. This choice is due to the fact 
that on the one hand the path of the OT light in the ISM is 
longer than that of field stars, and therefore the induced polar- 
ization may be different, and on the other hand the host galaxy 
ISM may also contribute to the induced polarization. Since the 
total induced polarization can be dealt with using a single set 
of q and u Stokes parameters, we prefer to use here the uncor- 
rected data (as previously done with GRB 021004, Lazzati et 
al. 2003), bearing in mind that the average Stokes parameters of 
field stars are qisM = 6.22 x 10~ 3 and uism = —3.95 x 10~ 4 
(Gorosabel et al. 2003). The second difference with respect to 
the dataset of Gorosabel et al. (2003) is that the observations 
have been binned in time in order to increase the signal-to- 
noise 1 . In particular, the first 4 data points (with 750 s exposure 

1 Polarization points that have been binned were consistent with 
each other. 



each) in Tab. 1 of Gorosabel et al. (2003) have been averaged 
into a single point, as well as the remaining 3 points of the first 
observing night (300 s exposure each). Our dataset is presented 
in Fig. [5] Even though we present polarization and position an- 
gle data, the fits were performed in the Stokes parameter space, 
where the uncertainties have a Gaussian distribution. 

Fit results are reported in Table ^ f° r the models de- 
scribed above. All the models have been fitted for evolution in 
a uniform ISM as well as for a wind environment. In the case 
of HHJ models, the possibility that the jet undergoes sideways 
expansion at the internal sound speed was considered as well. 
In addition to the jet parameters, the possible contribution of a 
polarizing ISM with free properties was considered. The result- 
ing fit was always consistent with the q\sm and uism derived 
from the field stars. We have therefore frozen the induced po- 
larization to the Galactic value in all the fits, in order to increase 
the number of degrees of freedom. This result is also consistent 
with the low upper-limit for reddening derived by Covino et al. 
(2003a) from U BVRIJHK quasi-simultaneous photometry 
and with the spectropolarimetry of Barth et al. (2003). In order 
to check also for the possible presence of an external ordered 
magnetic field (Granot & Konigl 2003) we have allowed the 
ISM Stokes parameters to become comparable to the measured 
polarization. This did not lead to a significant improvement of 
the fit. 

Fits with a hydrodynamic homogeneous jet yield always 
non-acceptable \ 2 values. This is mainly due to the fact that a 
minimum of polarization in coincidence with the break time is 
not observed, nor a 90° rotation of the position angle between 
the data taken across the jet break time. A reasonable (even 
though still not formally acceptable) fit can be obtained only in 
two limiting cases, which are not likely. If the jet break time is 
allowed to become smaller than 0.33 d, which is not consistent 
with the break time detected in the light-curve (Gorosabel et 
al. 2003) or if the line of sight is allowed to be slightly outside 
the jet edge in a non sideways expanding beam, the measured 
polarization points lie in a time interval where the angle rota- 
tion is not expected. In this case reasonably good fits can be 
obtained (with reduced % 2 values of ~ 2 3). The best HHJ fit 
is shown with a thin dashed line in Figs.|5]and[6] 

The fit with a structured hydrodynamic jet gives better re- 
sults. This is mainly due to the fact that in this case the position 
angle of polarization is constant, and the polarization curve has 
a maximum in coincidence with the break time. In this case a 
good fit is obtained (x 2 /d.o.f.=10.36/6) for a jet expanding in a 
uniform medium observed very near to its core, consistent with 
the small measured opening angle from the light-curve break 
time (Covino et al. 2003a; Rossi et al. 2002). In this case the 
break time can be measured also from the polarization points, 
with a result that is consistent with the limits from the light- 
curve. A moderately ordered magnetic field is also required. A 
fit with the same model expanding in a wind environment does 
not yield an acceptable fit. The best fit is shown with a thick 
solid line in Figs. |5] and |6] 

The best fit of the whole set is obtained for an MHJ jet, ir- 
respective of the environment structure. The ISM and Wind fits 
are shown with a thick dash-dot-dot-dot and long-dashed lines, 
respectively, in Figs. |3]and|6] Even though the fit is good, it is 
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Fig. 5. Polarization (upper panel) and position angle (lower 
panel) data for GRB 020813 (Gorosabel et al. 2003). Different 
curves refer to different models, as indicated in Fig. [5] Models 
yielding an acceptable fit are plotted with a thick line, while 
non-acceptable models are shown with a thin line. The gray- 
shaded area shows the acceptable range for the jet break time. 
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Fig. 6. Same as Fig.|5]but for the polarization only. The x and 
y axes have been expanded in order to show how the models 
differ at early times, where more data are needed in order to 
clearly understand the structure of the jet and of the magnetic 
field. The gray-shaded area shows the acceptable range for the 
jet break time (Gorosabel et al. 2003). 



Table 1. Fit results. 



Model (a) 


b (b) 


Ac) 
t iot 


#obs/#jct 


aid) 
u OT 


X 2 /d.o.f. 


HHJ ISM 


0.35 


0.33 


0.8 


151.0 


84.7/6 


HHJ Wind 


0.32 


0.33 


0.8 


151.0 


119.2/6 


HHJ ISM SE 


0.63 


0.33 


0.5 




194.0/6 


HHJ Wind SE 


0.57 


0.33 


0.5 


150.0 


194.2/6 


HSJ ISM 


0.475 ± 0.02 


0.29 ± 0.07 


2.0 ±0.1 


152.0 ± 1.2 


10.36/6 


HSJ Wind 


0.46 


0.33 


2.0 




18.7/6 


MHJ ISM 


< 0.6 (e) 


0.7 (e) 


0.4 (e) 


152.0 ± 1.2 


8.9/6 


MHJ Wind 


> 0.5 (e) 


0.38 (e) 


0.13 (e) 


152.0 ± 1.2 


8.7/6 


MSJ ISM 


0.185 


0.33 




152.0 


15.8/7 


MSJ Wind 


0.18 


0.33 




152.0 


29.3/7 



( a > Acronym of the considered model. Models yielding an acceptable fit have been highlighted in bold. Models with SE have sideways expansion 
included in the computation. 

Alignment of the magnetic field (see text). 
*- c ' The jet break time (in days) is constrained to be consistent with what measured from the light-curve (Gorosabel et al. 2003). In most cases 
the break time is forced to be the smallest possible (0.33 days). 
^ Position angle of the intrinsic OT polarization before the propagation in the ISM. 

' e ' Despite the goodness of the fit it was impossible to derive meaningful error intervals for these parameters (see text). 



not possible to independently constrain the parameters, since 
the covariances are strong and allow one to find always a min- 
imum of the x 2 within 1 from the absolute minimum. 

Finally, the MSJ model yields a reasonable fit, even though 
not formally acceptable. The fit is shown with a thin dash-dot 
line in Figs.|5]and|5] Taking into account the oversimplification 
with which this model has been computed, we do not consider 



the lack of a formally adequate \ 2 value a condition strong 
enough to reject the model. 

4. Summary and Discussion 

We have presented a comprehensive modelling of the polar- 
ization curve of GRB 020813. This burst is particularly suited 
for polarization studies since it has an extremely smooth light 
curve (Gorosabel et al. 2003; Laursen & Stanek 2003). This is 
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an important parameter, since any complexity in the light curve 
is likely associated with the breaking of the fireball symme- 
try, introducing a random fluctuation in polarization that can- 
not be predicted a priori by any model (Granot & Konigl 2003; 
Lazzati et al. 2003). The polarization curve of GRB 020813 
(Gorosabel et al. 2003) is one of the most extended published 
to date 2 , certainly the most complete associated to a smooth 
light curve, and is characterised by a constant position angle 
and a smoothly decreasing degree of polarization. Importantly, 
the data encompass the break time of the light-curve, which is 
a critical time where different models make markedly different 
predictions (Rossi et al. 2004). 

We first discuss the magnetic patch model (Gruzinov & 
Waxman 1999), in which polarization is due to the non-perfect 
cancellation of highly polarized radiation coming from a large 
number of ordered magnetic field domains independent from 
each other. This model predicts a strong flickering of the posi- 
tion angle which is not observed and can therefore be rejected 
based on our dataset (see also Barth et al. 2003 and, for the case 
of GRB 030329, Greiner et al. 2003). 

We then model the polarization curve according to the pre- 
dictions of several different models. We find that a homo- 
geneous jet in which the magnetic field is shock generated 
(Ghisellini & Lazzati 1999; Sari 1999) cannot fit the data, since 
neither a minimum of polarization nor a rotation of the position 
angle are present in the data. On the other hand a structured 
model, in which the core of the jet is more energetic than its 
wings, can successfully reproduce the data, and predicts a jet 
break-time in agreement with what measured from the light- 
curve. We also compute models for a magnetised jet, in which 
the magnetic field has a non negligible toroidal component. We 
find that the best fit is obtained by a homogeneous magnetised 
jet, even though the dataset is not extensive enough to mean- 
ingfully constrain its parameters. We also consider a force-free 
magnetic bubble (Lyutikov et al. 2003), which is analogous to a 
structured jet but for the presence of a toroidal magnetic field. 
The fit in this case is not formally successful, but we cannot 
rule out the model since the details of how the toroidal mag- 
netic component mixes with the shock generated one have not 
been deeply investigated. 

All the considered fits require a non fully ordered magnetic 
field, even though the degree of order is always substantial 
(> 40%). This is in disagreement with what is found in the 
prompt emission of GRB 021206 where a fully aligned mag- 
netic field has to be considered in order to reproduce the obser- 
vational constraints 3 (Coburn & Boggs 2003; Lyutikov et al. 
2003; Granot 2003; but see also Nakar et al. 2003 and Lazzati 
et al. 2004). This indicates that the magnetic field geometry 
is different in the two epochs and seems to confirm the idea 
that the material responsible for the prompt emission is not the 



2 A more extensive polarization covering has been performed on 
GRB 030329 (Greiner et al. 2003). This burst, however, has a complex 
light curve which prevents the comparison of the polarization with 
models, as for the case of GRB 021004 (Lazzati et al. 2003). 

3 Note however that the result has been heavily criticised on data 
analysis grounds (Rutledge & Fox 2004) and any implication should 
therefore be taken with caution. 



same one that produces afterglow photons. Our best fit model 
requires however a mixing of the two, or at least of their elec- 
tromagnetic components. 

Unfortunately the quality of the data is not good enough 
to allow us to constrain the models and/or the assumption on 
which they are based any further. Robustly and independently 
of the model assumptions we can conclude that: i) the magnetic 
field responsible for afterglow emission is not simply random 
since a net polarization signal is produced; ii) the difference 
between the magnetic field components parallel and perpen- 
dicular to the shock front is small but non negligible, but it is 
not possible to understand which of the two dominates over 
the other; iii) a standard homogeneous jet with shock gener- 
ated field can be ruled out, since no angle rotation is detected 
between observations before and after the jet break 4 . On the 
theoretical side, we show that the early time behaviour of po- 
larization is crucial to understand the magnetic field configura- 
tion: a large scale field will produce large polarization at early 
time, while a small-scale field will produce polarization only 
within a couple of decades in time from t ~ tj et . 

We conclude by commenting on how future observations 
may help to gain a deeper insight on the structure of GRB jets 
and their magnetic field. As shown in Fig.|6j the various models 
differ substantially at early times. In particular, magnetised jets 
predict high polarization at early times, while unmagnetised 
models predict null polarization. Early polarization measure- 
ments of afterglows with smooth light curves will therefore be 
fundamental to pin down these important jet parameters. 

We shall also conclude with a word of caution. Dense sam- 
pling of the polarization curve of GRB 030329 (Greiner et al. 
2003) has revealed that the polarization curve associated to a 
complex light curve afterglow can be much more complex than 
what predicted by any of the models we discussed here, and 
that fluctuations in the light and polarization curves may not 
be strictly correlated. In the case of GRB 020813 the light and 
polarization curve sampling is not as dense, and therefore even 
if all data are consistent with a smooth evolution, short time- 
scale variability in both the light and polarization curve may 
have been missed 5 . We therefore strongly recommend an ade- 
quate sampling of GRB afterglow light curves in polarimetric 
mode in order to allow for a more robust comparison of the data 
with models. 
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A constant position angle does not imply necessarily an inho- 
mogeneous jet with the structure adopted in this paper. Any jet 
with a brighter core and dimmer wings will produce qualitatively 
similar polarization curves. 

5 However Lazzati (2004) showed that a degree of fluctuation as 
large as that in GRB 030329 would have been easily detected in 
GRB 020813. 
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